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Abstract 

The linking and writhing numbers are key quantities when characterizing the struc- 
ture of a piece of supercoiled DNA. Defined as double integrals over the shape of the 
double-helix, these numbers are not always straightforward to compute, though a sim- 
plified formula exists |13] . We examine the range of applicability of this widely-used 
simplified formula, and show that it cannot be employed for plectonemic DNA. We 
show that inapplicability is due to a hypothesis of Fuller theorem [13] that is not met. 
The hypothesis seems to have been overlooked in many works. 

PACS numbers: 87.10.-e, 87.14.gk, 02.40.-k, 87.15.-v 



1 Introduction 

The double-helix structure of DNA is the source of many complications in its in-vivo func- 
tioning during condensation/decondensation, replication, or transcription. For example the 
necessary unzipping of the molecule during transcription induces torsional strain along DNA. 
The development of magnetic tweezers [311 135] or the use of quartz cylinders in optical tweez- 
ers [1] allow researchers to investigate the in-vitro response of DNA molecules to torsional 
stress. Studies of the behaviour of this twist storing polymer are not just a game for physi- 
cists as it has been clearly established that, e.g. the assembly of RecA could be stalled 
by torsional constraints [41], or that the rate of formation and the stability of the complex 
formed by promoter DNA and RNA polymerase depends on the torque present in the DNA 
molecule ^B]. On the theoretical side, matters are made difficult by the nonlocality of the 
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topological property that is associated with the torsional constraint: the link. The two sugar- 
phosphate backbones of DNA have opposite orientation and the ends of a double-stranded 
DNA molecule can only be chemically bound in such a way that each strand joins itself, 
thereby yielding two intcrwound closed curves (no Mobius-like configuration can exist). For 
a circularly closed DNA molecule, the link is the number of times one of the sugar-phosphate 
backbone winds around the other. Once the three dimensional shape of the molecule is pro- 
jected onto a plane, the link is given by half the number of signed crossings between the 
two backbones. This quantity is best seen as the number of turns put in an initially planar 
plasmid (a piece of circularly closed double-stranded DNA) before closing it. Link has been 
shown to consist of two parts: the two sugar-phosphate backbones of a plasmid can be linked 
because (z) the plasmid lies in a plane but the base-pairs are twisted around their centre line 
(the curve joining the centroids of the base-pairs) and/or (ii) the centre line itself follows 
a writhed path in space. In general the two possibilities coexist and the link Lk of a DNA 
molecule is the addition of the two quantities [71 HH E] : 

Lk^Tw + Wr. (1) 

The twist Tw is a local quantity in the sense that it can be computed as the single integral 
of the twist rate t(s): 2t: Tw = T{s)ds, where s is the arclength along the centre line 
and L the total contour length of the molecule. In elastic DNA models, the twist rate r(s) 
is normally coupled to mechanical quantities (e.g. the torque) characterizing the molecule. 
Contrary to the twist, Wr is a global property of the centre line of the molecule. We first 
define the directional writhe of a closed curve F. Consider a closed curve in 3D and project 
this curve, along a certain direction, on a plane. The number of signed crossings seen in the 
plane is the directional writhe for that direction. One could then reiterate the procedure with 
different directions of projection and compute the directional writhe for each direction. The 
average value obtained for the directional writhes, when all directions are considered, is the 
writhe Wr{T) of the 3D curve F [H]. 



2 Global and local writhe formulas 

The following double integral has been introduced by Calugareanu [H [7] and White [42] : 



where r{s) is the position of points on F and t{s) is the unit tangent to F. As soon as a 3D 
curve does not self-intersect, the writhe Wr{T) of the curve is given by Eq. ([5]): Wr{T) = 
Wr^^ (T) if r(s) ^ ^(s') Vs, s' with s ^ s' [^. The computation of the double integral of 
Eq. ([2) is analytically hard and numerically time consuming to perform. An important 
result though enables one to reduce the double integral to a single integral, provided several 



^The double integral may be shown to converge for self-crossing curves provided the tangents at the 
intersection points are not aligned 1331 . 
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hypotheses are fulfilled [T^ . It is the main purpose of this paper to show that these hypotheses 
are not met in the case of plectonemic DNA. 

Fuller's theorem [13] states that the writhe of the curve F can be computed by considering 
the writhe of a reference curve Tq (which should be known or easy to compute) and the 
continuous deformation (a homotopy) morphing Fq to F: 

WriT) ^Wr^{T,To) Wr{To) + ^FuiT,To) (3) 

where Fu(r,Fo) = P J^^^^ . A + , 

Jo I + to{(7) ■ t{a) da 

where to[(j) =^ {dro{(T)/da)/\dro{(T)/d<T\ is the unit tangent to Fq and where a G [0,(7^] is 
a common parametrization for both curves, though not necessarily the arclength along F 
or Fq. The total contour length of F is L. The curve Fq is normally chosen to be a fairly 
simple curve (e.g. planar with zero writhe) so the computation of Wr{T) boils down to 
estimating the single integral in Eq. Before doing so, one should verify whether the 

hypotheses of Fuller's theorem are fulfilled (this is the sense of the question mark in Eq.(I2])). 
The continuous deformation from Fq to F introduces a familly of curves F^ with A G [0, 1] and 
Fi = F. The unit tangent to F^ is t\{a) '= {dr\{a)/da)/\dr\{a)/da\. The first hypothesis 
is that none of the curves F^ self-intersects. The second hypothesis is that, and this is the 
point we want to emphasize, there should be no point along any of the curves F^ where 
t\{a) ■ to(cr) = —1. For each value of A, the (unit) tangent t\{a) with a G [0,(7^] defines 
a curve, called the tangent indicatrix, on the unit sphere. If for a certain A and a certain 
a we have t^^a) ■ tQ{a) = — 1 then, on the unit sphere, the point corresponding to t^i^) is 
antipodal to the point corresponding to to(5'). For brevity we shall also call antipodal the 
point r;^(a) or the curve itself. 

In models where self-intersection is prevented (e.g. models using hard- or soft-wall poten- 
tials of the molecule on itself) a continuous deformation free of self- intersection (an isotopy) 
may be easy to devise and the first hypothesis would be verified. However nothing in these 
models ensure that the second hypothesis (that no antipodal points exist for all A G [0, 1] 
and all a G [0, cl]) is met and consequently Fuller's formula ^ cannot be used to compute 
the writhe. We nevertheless remark that even in the case of self-intersection or antipodal 
points Fuller's formula is always correct modulo 2: the integer, and most important, part 
of Wr(r) will not be correct, but the fractional part will. Each antipodal point present in 
the continuous deformation from Tq to F introduces a shift of two units between the actual 
writhe Wr(r) of the curve F and the value W^r^(F, Fq) 1 . Consequently, if there are m 
antipodal points along the deformation A G [0, 1], the computed value Wr^ could be as far 
as 2m from the correct value Wr, but the fractional part will be accurate. 

This discrepancy between Wr^ and the actual writhe Wr has already been pinpointed 
in the case where DNA is treated as a fiuctuating chain under low twist (i.e. without plec- 
tonemes) [28l ESI HO]. Nevertheless we have found a certain number of references where 
Fuller's formula is used [IHl [III S l [3 [HI [11 III [II |1 111 [13 III |13 113 113 |M IS In 
most cases the hypotheses of Fuller's theorem were not checked, if only mentioned. In some 
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works the formula is used in a scheme that provides an estimate for the torsional stiffness 
of the DNA molecule [H [3] and this has been shown to lead to incorrect results [27l [26l EH] • 
Nevertheless, in some few papers [301 HOI Ell tUji the formula was used to asses the writhe 
of DNA configurations under high stretching force and low twist, in which cases antipodal 
points are absent and the formula is correct. 

3 The writhe and Unk of DNA molecules in magnetic 
tweezer experiments 

During a force-extension experiment on a single DNA molecule thermal agitation deforms the 
molecule whose shape locally adopts random directions (Fig. [H left). In the absence of (or 
under low) twist, DNA is modeled as a worm-like chain [T^]; the path followed by its centre 
line in space looks like a (directed) random walk [20]. In such configurations, the writhe is 



magnetic 




Figure 1: The shape of a DNA molecule subjected to a vertical pulling force follows a 3D 
curve which is locally made random by thermal agitation. As explain in [28j , when the force 
is small or moderate, points where the tangent to the molecule is oriented downward can 
appear, e.g. point A. These points on the curve are antipodal to corresponding points on the 
reference curve. 

usually evaluated using Fuller's formula ([3]) with the reference curve shown in Fig. [H left. 
Since the writhe is classically defined for a closed curve, the reference and actual curves are 
all closed by imaginary C-shaped curves (dashed in Fig. [T]) that connect the top to the base 
of the configurations. (Other choices of closures and interferences due to the closure in the 
computation of the writhe are discussed in [28l [33l EH Ell HO] . ) Strictly speaking the curve F 
for which we compute the writhe consists of two parts: the closure and the part corresponding 
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to the molecule. As the closure remains unchanged in the continuous deformation A G [0, 1] 
its contribution to Fuller's integral is zero. Consequently we will forget about the closure 
and call F the (open) curve described by the molecule only. Fhe reference curve Fq is then 
the z axis, and Fuller's formula ^ becomes: 



T4^r^(F,Fo) = T4^r(Fo) + 

ZTT 



1 



• f (s) as 



= 



1 

2^ 



(1 - cos a{s)) — as , 



(4) 



if the tangent t(s) = (sin 61 cos V', sin sin ^, cos 0) is parametrized by Euler angles 0(s) G [0,7r], 
and 'ijj{s) G [0, 27r). As stated above, Eq. (j4|) is only valid if the curve F can be deformed 
into the reference curve Fq without passing through configurations Fa having their tangent 
vector t\{a) facing -e^: t\{ij) ^ — for all A, a. We remark that in the present case where 
Fq is the z axis an antipodal point is a point where the curve t\ (cr) passes through the south 
pole of the unit sphere: 0{a) = tt. 



c 



4 





Figure 2: Configurations A, B, C, D, and E (see Fig. |4|) obtained by numerical continuation 
[23] . Configurations B and D each have an antipodal point to the reference curve of Fig. [TJ 
These antipodal points are located at the middle point of the end-loop of the plectonemic 
structure. 



In the case where the actual configuration F exhibits a single or multiple antipodal points 
it has been proposed in [3] to evaluate the writhe from Eq. (j4]) with Euler angles defined 
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on a truncated unit sphere: 6 would not be allowed to reach tt and hence antipodal points 
would be avoided (the curve T so defined would be very near the real curve F, hence the 
writhes would almost be the same). We stress that this process is not sufficient as even for 
a given curve with no antipodal points, it is not clear whether Eq. ^ is valid or not. In 
order for Eq. ([¥]) to be valid, one has to exhibit a continuous deformation from the z axis to 
the curve F that is entirely free of antipodal points. There are many cases where the actual 
curve F is free of antipodal points, but where a deformation free of antipodal points does not 
exist. Consequently in these cases Eq. @ yields an incorrect result (unless antipodal points 
of opposite signs cancel out). Nevertheless, it has been noted that under high stretching 
force (e.g. F ^ h pN) and low torque the DNA molecule is almost straight and no such 
antipodal points exist [22]. In such a case Fuller's formula is correct provided no plectonemes 
are present, as we will see now. 





A 


B 


C 


D 


E 


Wr 


0.70 


0.99 


1.09 


2.92 


3.45 




0.70 




-0.91 




-0.55 


Lk 


4.46 


1.84 


4.15 


4.23 


4.75 




4.46 




2.15 




0.75 



Table 1: Writhe (first line) and Link (third line) of the configurations of Fig. [21 computed 
by continuation or use of the double integral (Eq. ([2|)). The second line is computed from 
Eq. ([4|) and the last line is computed from Eq. ([6]). Formulas ([4|) and ([6]) are not applicable 
on configurations B and D that each have an antipodal point, and yield incorrect results for 
configurations C and E. 



In magnetic tweezer experiments, when a large amount of turns are put in (by rotation 
of the magnetic bead around the z axis) , the DNA molecule reacts by forming plectonemes 
[36j . The number of turns n imposed on the magnetic bead is given by the Link Lk of the 
molecule. We now show that the presence of plectonemes in the supercoiled configuration 
prevents the existence of a deformation, from Fq (i.e. the z axis) to F, that is free of antipodal 
points, and consequently forbid the use of Eq. Q . 

In terms of the Euler angles the twist of the molecule can be computed as 

+ -"?)*. P) 



27r Jq \ds ds ^ 

where 4>{s) is the third Euler angle (see e.g. [39]). Using Eqs. ([4|) and ([5]) we obtain 

1 /-^ d 



Lfc^^-y^ -(0-^v.)rf.. (6) 

This formula usually is the starting point for computations of the link of supercoiled configu- 
rations, see e.g. [3] or [TU]. We show here that it yields incorrect results when plectonemes are 
present. We have performed computations to model the elastic response of a twist storing 
filament subject to tensile and torsional constraints and we quantitatively reproduced the 
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Figure 3: Tangent indicatrices of configurations A, B, C, D, and E of Fig. ^ An antipodal 
point arises when the curve passes through the south pole of the sphere, i.e. for configurations 
B and D. 



plectonemic regime characterised by the linear decrease of the end-to-end distance of the fila- 
ment as a function of the number n of turns put in [23| . The plectonemic configurations were 
computed numerically using a continuation algorithm and n — Lk was computed by conti- 
nuity so that no integer number of turns is missed (we also performed numeric integration of 
the double integral of Eq. ^ for the writhe of the configurations, with closures, and always 
obtained consistent results) . These shapes serve as an illustration for the computation of the 
writhe and consequently are used for their geometry only. The fact that they are mechanical 
equilibria is not relevant here. The continuous dark curve in Fig. H] shows an output of the 
numerics, drawn in the {Lk, Z) plane, where Z is the vertical extension of the molecule. The 
curve starts at point O which corresponds to a straight and twisted configuration. The path 
is then monotonically decreasing in Z. We have selected five configurations A, B, C, D, 
and E which are drawn in Fig. [H together with their corresponding tangent indicatrices (see 
definition in Section [2]) drawn in Fig. [31 Configurations B and D each comprise an antipodal 
point located at the middle point of the end-loop of the plectonemic structure. This can be 
verified in Fig. \^B and D where the tangent indicatrices pass through the south pole of the 
unit sphere. Using the same geometric configurations, we compute Lk^ from Eq. ([SJ and 
we plot the corresponding curve, in gray, on the same diagram. We see that each antipodal 
event introduces a shift of two units in the gray curve, which is consequently broken. As 
expected this confirms that Eq. ^ is only valid modulo 2 and hence should not be used 
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Figure 4: The dark continuous curve is obtained from numerical continuation [23] and the 
broken gray curve is obtained from Eq. Five configurations A, B, C, D, and E are 
selected and plotted in Fig. [2] together with their tangent indicatrices in Fig. [31 The inset 
shows a zoom around point B, the first antipodal point, where the gray curve jumps by an 
amount of -2 units. From O to B the two curves coincide. 
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to estimate the link of plectonemic configurations. The writhe and link of the five example 
configurations are given in Table [T] where they are compared with Wr^ and Lk^ given by 
Eqs. (jll) and For configuration E, which is separated from the reference curve by two 
antipodal events, we see that Wr^ (resp. Lk^) is 4 units away from the correct Wr (resp. 
Lk) value. This configuration E is an illustration of the fact that Fuller's formulas (Eq. Q 
or ([6])) can be wrong even for a configuration that does not comprise any antipodal point, 
which is clearly illustrated in Fig. ^E where we see that the tangent indicatrix is nowhere 
near the south pole. 

We first comment on the size of the gap between the two curves. A plectonemic DNA 
configuration F that shows (positive) crossings on a lateral projection has Wr{T) ~ +Nx 
[55] . If we continuously deform this configuration to the reference curve Fq of Fig. [T] left by 
unwinding the plectonemic region, an antipodal point arises each time the tangent at the 
apex of the terminal loop (at the end of the plectonemic region) is facing downward. This 
happen Nx/2 times. As proved in [Tj the presence of m = Nx/2 antipodal points leads to a 
discrepancy in Fuller's formula of 2m: \Wr — M^r^j < 2m. Since 2m = ~ Wr, we have 
that: < Wr^ < 2Wr, which corresponds to an error of up to 100%. This is apparent in 
Fig. [4] where the broken gray curve stays near the vertical axis while the (continuous) dark 
curve monotonically increases in Link. 

Second we note that the discrepancy between the two curves occurs shortly after self- 
contact has started in the filament, for n ~ f , see inset of Fig. U) Plectonemic structure may 
exist for small number of turns n, provided the pulling force is not too large. On the other 
hand, a large pulling force does not rule out the occurence of plectonemes, provided that n 
is large enough. This leaves a small parameter regime (large pulling force, low number of 
turns) where plectonemes are absent. When the two sources of discrepancies are considered 
(random walk antipodal points ^Hl HE] , and plectonemic antipodal points) one sees that the 
use of Eq. ^ (resp. Eq. ©) to compute the writhe (resp. the link) in a model for DNA 
under tensile and/or torsional stress is to be avoided unless n ~ and the tensile force is 
large. 

4 Discussion and concluding remarks 

We summarize here few properties of the quantities VFr'^'^(F), VFr^(r,Fo) which do not 
always give the writhe of a curve F. The quantity Wr^^ (T) is a function of the curve F 
only, whereas the quantity Wr^{T, Fq) also depends on the reference curve Fq. The quantity 
IYj-CW yigi(jg tiig correct value for the writhe Wr of a closed curve as soon as the curve is 
not self-intersecting. Along a continuous deformation F^ with A G [0,1] (with Fi = F), the 
quantity Wr^^iTx) jumps by two units when the curve F^ intersects itself (say at A = 1/2). 
Since the writhe Wr also jumps by two units, the quantity Wr^^ is equal to Wr before and 
after the self-crossing event, i.e. Wr'^^iTx) = Wr{Tx), VA 7^ 1/2. Now the quantity Wr^ 
has no such discontinuity: along a continuous deformation where the curve self-intersects, the 
writhe Wr will jump by two units, but the quantity Wr^ will stay continuous. This means 
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that the quantity Wr^ no longer yields the correct value for the writhe after the self-crossing 
event, Wr^{rx,To) = T^r(rA), VA < 1/2 but Wr^{rx,ro) ^ Wr{Tx), VA > 1/2. The same 
is true for antipodal points. Along a continuous deformation Tx with A G [0, 1], the quantity 
Wr^ [Fx, To) jumps by two units when the curve Fx has an antipodal point (say at A = 1/2) 
with regard to the reference curve Tq. On the other hand the writhe Wr stays continuous. 
This means that the quantity Wr^ no longer yields the correct value for the writhe as soon as 
an antipodal event happens, Wr^{Tx,To) = Wr{Tx), VA < 1/2 but Wr^ {rx,To) ^ T^r(rA), 
VA > 1/2. This is an important point and many authors seem to believe that the quantity 
W^^'^lrAjTo) only has problems for configurations actually comprising an antipodal point 
(A = 1/2 in the above example). Therefore the quantity M^r^(r,ro) may not be equal to 
the writhe H^r(r) even for curves that do not comprise any antipodal point. This make the 
use of l¥r^(r,ro) uneasy, as one has to first verify the absence of antipodal events in the 
entire continuous deformation Fq — * F . On the contrary the use of Wr'~^^ (F) is much easier 
in the sense that one just has to check that the actual curve F does not self-intersect. In this 
sense the quantities VFr'^^(F) and M^r^(F,Fo) do not suffer from the same pathologies in 
the computation of the writhe, contrarily to what is claimed in [29 . Another consequence 
is that, when sampling DNA configurations to construct a statistical ensemble and compute 
writhe averages and fluctuations, it is not enough to introduce, as was done in [5, a small 
forbidden region around the south pole of the unit sphere to ensure that Wr^ i^, Fq) yields 
a correct value. In fact many of these sampled configurations, even free of antipodal points, 
are configurations that suffer the same problems as the configurations with A > 1/2 above: 
iyr^(F,Fo) 7^ M^r(F), as numerically verified in [^[^. 

Finally we want to point out the following property. We saw that if a continuous de- 
formation Fa with A G [0,1] contains an antipodal event, then M^r^(F,Fo) does not yield 
the correct result after the antipodal event. Now one can argue that yet another continuous 
deformation, free of antipodal events and self-crossings, may exist (with the same initial and 
final curves) and that in this case Wr^ {T^Yq) would yield the correct result. This is not the 
case since, as we show in Appendix B, as soon as a corrupted deformation exists in between 
two curves Fq and F, then all continuous deformations are corrupted and 14^r^(F,Fo) defi- 
nitely yields an incorrect result. This means that the reference curve Fq cannot be used. One 
is bound to find another reference curve or to use the double integral Wr'^^ . It appears that 
an easy way to verify the applicability of Fuller formula Q between an actual curve and a 
reference curve is to consider any convenient deformation that avoids antipodal events. This 
can always be done if one allows for self-crossings. Then formula ^ is applicable if and only 
if the net sum of self-crossings is zero. 

In the case of numerical computations, either dealing with continuation of mechanical 
equilibria or statistical ensembles of configurations, the writhe can be assessed in an accurate 
way by using both the double integral ^ and Fuller integral (H]) in a cooperative way. The 
double integral being time consuming to evaluate, one can in a first step discretize it |37] 
with a reduced number of elements so that it produces an approximate result that only has 
to be accurate up to ±1 (one still has to estimate how many elements are needed to obtain 
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such an accuracy [4]). Then in a second step, FuUer integral is used to refine the result. 
Even if Fuller integral may be off by several integers, its fractional part is correct. In this 
scheme, the double integral yields the integer part of the writhe and Fuller integral yields 
the fractional part of the writhe. Moreover in cases where one knows the correct value of 
the writhe of a nearby configuration (e.g. a predecessor configuration, one move away in a 
Monte-Carlo scheme) this value can be used as the approximate result of the first step and 
one only has to compute Fuller integral to obtain a correct and accurate value of the writhe. 
(The usual assumption that there is no self-crossing between the predecessor and the actual 
configuration still holds.) 

In conclusion we have shown, by producing counter-examples and explaining the under- 
lying causes, that formulas ^ and ^ cannot be used to compute the writhe and the link 
of supercoiled DNA configurations encountered in magnetic tweezer experiments. 
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6 Appendix A : Fuller's formula in a detailed example 

In this appendix we compare, for a given example curve, the value computed from the 
double integral Wr^^ to the value computed from the single integral Wr^ . We start by 
introducing the parametrization of a circularly closed supercoiled plasmid (see Fig. [5]). The 
parametrization is divided in four parts, called A, B, C, and D. Part A is a right-handed 
helix of radius unity, total arclength La, and helical angle 9a (with < Oa < 7r/2): 

(simpA{sA) \ I sin e A cos iI:a{s a) \ 

-cosiI}a{sa) and <A = — = sin 9 a s\n if: a{s a) (7) 
S A cos 9 A j V cose A ) 

with SA G [0, La\ and V'A(syi) — sa sin 6*^. We note AV'a =^ La sin^^. We have |tA(sA)| = 1 
for all Sa, and consequently sa is the arc length in part A. The end loop B connects the two 
helices A and C: 

(sin cossb -I- cos A^/;^ sin^yi sinsB cos s_B \ 
— cos A'f/'^ cos SB -I- sin A-0^ sin 6*^1 sin SB cos SB and — (8) 

LACOseA + {s\n9A-sin9B{sB)) /{1 + 29a/tt) ) '^^^ 

with SB e [0,7r], and 9b{sb) = 9a — sb[1 + 20^/7r). For this part B, \tB{sB)\ i= 1- 
Part C is a right-handed helix of radius unity, total arclength Lc — La and helical angle 
9c — — TT — 9a = ''^ — 9a mod2'K: 

I -sinV'c(sc) \ dr ( sin 6*0 cos ?/'c(sc) \ 

rc = I cosi)c{sc) I ^'^'^ '^'^ ^ ^ I sin 6'c sin '0c(sc) I (9) 

y La cos 9a + sc cos 9c j \ *^°s 9c j 
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with sc G [0, Lc] and i^dsc) — ^i'A — sc shiOc- We have \tc{sc)\ — 1 for aU sc, and 
consequently sc is the arclength in part C . The end loop D closes the curve: 

(sin 9 a sin s d cos so \ 
cossd and <r) = — (10) 

(sin0z3(sz5)-sin0A)/(l + 2^?A/7r) / 

with Sd S [0, tt], and 9d{sd) = — tt — 0a + sd (1 + "^Qa/t^)- For this part _D, |tr)(sD)| 7^ 1. 
We call r the union of the four parts: V — Ta U Tb U Tc U F/). The curve has continuous 
derivatives, and finite jumps in its curvature. 




Figure 5: Circularly closed supercoiled plasmid corresponding to curve F — FaUFbUFcUFd 
with 0A = 7r/6 and La = Stt. 



Analytical calculation 

The evaluation of the double integral of Eq. ^ in the limit La > 1 yields [Ml [551 [M] : 

M^rC^^(r)^^°°-^^sin0ACOS0A. (11) 

TT 

The minus sign is due to the fact that the helices A and C are right-handed. 

We then compute ^^^^(r.ro), given by Eq. ([3]), using the reference curve Fq of Fig. El 
top. The continuous deformation between Fq and F is then parametrized by (?Aa £ [0, 6'a]. 
The reference curve of Fig. [5] top comes to mind naturally and makes calculations easiest, 
but yields antipodal events and hence an incorrect result, as we shall see now. The writhe of 
the reference curve Fq is zero. Consequently we only have to compute the integral Fu{T^ Fq) 
of Eq. ^ over the four different parts A, B, C, and D. Clearly when the plasmid is long 
enough, i.e. when La is large, the contribution to Fu{T, Fq) of the helical parts A and C 
become dominant. Indeed their contribution scales with La while the contribution of the 
end loops B and D remains bounded. Consequently we focus on the contributions of the 
helical parts A and C, for large values of La, that is we look at the limit Fu(F,Fo) — > 
Fua(F,Fo) + Fmc(F,Fo) when La 00 5J. For the helical part A, the corresponding 
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tangent of the reference curve is t^o = (0,0, 1). Fuller integral for part A is then: 

FuA ^ / — — (1 — cos0A)ds^ sin^^ (1 — cos^yi) . (12) 
Jo dsA 

For the helical part C, the corresponding tangent of the reference curve is tco = (0, 0,-1). 
The Fuller integral for part C is then: 

Fuc= — P- {1 + cos 0c)dsc =smecil + cos Be) Lc ^ FuA (13) 
Jo dsc 

Neglecting the end loop contributions we arrive at: 

w^r^(r,ro) = w^r(ro) + ^Fu(r,ro) (14) 

= ^{Fua + Fub+Fuc + Fud) (15) 
zvr 

~ ^{FuA + Fuc)^— sine a{1- cos 9 a) (16) 
zvr TT 

We first remark that W^r^(r,ro) > ahhough Wr^^ (T) < (see Eq. (dU). The discrep- 
ancy Wr^(r) ^ Wr'-''^ (r) is due to the fact that in the continuous deformation 6a^ G [0, 9a] 
from the curve Fq of Fig. E] top to the plasmid of Fig. there are antipodal events. These 
events happen each time the middle point of the end loop B is pointing towards — Cj, and 
hence become antipodal to the corresponding tangent in Bq. As 9a)^ is increased, such 
an event happens once each time ^ipAx La sin^^^ increases by an amount of 2tt. In 
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the entire deformation there will be A^^/(27r) antipodal events. Each event introduces 
a shift between Wr^^iT) and W^r^(r,ro) of 2, the difference PFr^(r,ro) - Wr^^{T) 
should then be equal to 2 A'i/;yi/(27r) = (L^/tt) sin6'A, which can be verified by compar- 
ing Eq. (fTTj) to Eq. (fTH|) . A natural way to try to avoid the presence of these antipo- 
dal points is to rotate the reference curve, e.g. choose Tq' of Fig. [S] bottom as reference 
curve. In this case Fua = — cos^^sin^^ J^^^ smil^/{l + sin 0^ sin f/)) dip and Fuc = 
cos Oa sin 0A J^^"^ sin ip/{l — sin 9a sin ijj) dip ^ Fua ■ Here again the sum Fua + Fuc dom- 
inates Fu{T, Fq). Integration shows that Fua + Fuc — ^La sinf?^(l — cos 6 a) + vi^A, Ga), 
with < 1 — cos^A, which is not the correct result. Antipodal events are in fact 

still present with this rotated reference curve Fq and trying another rotation will not help, 
as explained in Appendix B. 

Numerical verification 

In Fig. [HI a right-handed circularly closed plasmid is drawn with 9a — 7r/6, La — 57r. The 
viewpoint selected to draw the curve in Fig. [5] is such that three (negative) crossings appear. 
For some other viewpoints, only two crossings appear. The writhe being the average number 
of signed crossings one sees from all possible viewpoints, it is relatively easy to convince 
oneself that its value for the curve of Fig. Olies in between -2 and -3. A numeric discretization 
scheme ^ of the double integral of Eq. ^ with 600 points yields Wr'^^{T) -2.2049. 
A numerical integration of Fu(F,Fo) = Fua + Fub + Fuc + Fud yields Wr^ {T,Tq) (= 
Wr^{T,T'o)) ^ -0.2050. The difference between Wr^^iT) and VFr^(F,Fo) is -2, which 
indicates that the net number of signed antipodal events is -1, as shown in Appendix B. 

7 Appendix B : The closed circuit theorem 

Antipodal points and self-crossings 

In Fig. [7] we show a family of curves along a continuous deformation. The deformation starts 
and ends with the same curve: we have a closed circuit. The 'stadium' shaped curve F^^o 
is continuously deformed through curves F^, as A increases from A = to A = 1, to the 
supercoiled plasmid Fi. The plasmid is then deformed back to the 'stadium' shaped curve, 
through another path, as A increases from A = 1 to A = 2. The values A = and A = 2 
refer to the same curve, the reference curve, which plays a special role. This reference curve 
is chosen to be fairly simple, in particular T/Fr(Fo) = 0. All the curves in the continuous 
deformation share a parametrization with parameter s (not necessarily the arclength) taking 
values from s = to s = L. We recall that an antipodal point is a point on a curve F^ such 
that the tangent t\{s*) to F^ at s = s* is aligned and in opposite direction with the tangent 
to{s*) to the reference curve Fq, i.e. an antipodal point is such that tx{s*) ■ to{s*) = — 1. 
An antipodal point is always defined with regard to a certain reference curve. Along a 
continuous deformation there may be A values where the curves T\ have antipodal points to 
Fq, and there may be A values where the curves Fa self-intersect. The closed circuit theorem 
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claims that in any closed circuit the (signed) number of self-intersections is equal to the 
(signed) number of antipodal points. For example, in the closed circuit of Fig. [7] there is 
one antipodal point and one self-crossing. This theorem was proposed and proven in |29] 
for a particular case of reference curve, a straight line. Here we give an outline of a proof 
which is valid for more general reference curves: Formula ^ has a ±2 discontinuity at each 
self-crossing event encountered in the closed circuit, but no discontinuity during antipodal 
events. On the contrary formula ^ has no discontinuity during self-crossing events but a 
±2 discontinuity at each antipodal event [l]. Formula ^ and ^ by definition agree when 
applied on the reference curve, and furthermore they may only differ by an even integer: 
Wr'^^ = Wr^mod 2. If we now follow the closed circuit, each formula will have its own 
discontinuities but both must eventually agree again at the end of the circuit. It follows that 
the number of discontinuities (respectively events) of one kind must be equal to the number 
of discontinuities (respectively events) of the other kind. 

No parallel between Wr^^{T) and W^r^(r,ro) 

The quantity Wr'~'^ (T) is not defined on curves that are self- intersecting, e.g. curve C in 
Fig. [3 The quantity Wr''{T,To) is not defined on curves F that have an antipodal point 
to Fp, e.g. curve A in Fig. \7\ These two properties, once considered with the fact that 
in any closed circuit the number of configurations with self-crossing is equal to the number 
of configurations with antipodal points, could lead one to believe that some parallel exists 
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between Wr^^ (T) and Wr^ (T.To)- In fact in [29 it is proposed that, in the computation 
of the mean writhe (or hnk) of a statistical ensemble of curves, avoiding self-crossings was 
equivalent to avoiding antipodal points. Moreover it was infered that for such statistical 
ensemble of curves, Fuller's formula, Eq. could be used safely provided that curves with 
antipodal point(s) were discarded when generating the statistical ensemble. We argue that 
this is not the case as Fuller's formula does not yield the right result for the (many) curves 
that have no antipodal point: e.g. curve B in Fig. [T] Again, Fuller's formula can only be 
used for a curve T^^ if one can devise a continuous deformation from Tj^ to Tq where none of 
the curves Fx with < A < A has antipodal points. For curve B in Fig. [7] this is not the case 
since the curve A has an antipodal point. Accordingly Fuller's formula does not yield the 
correct value for the writhe of curve B, as explicitly computed in Appendix A. In conclusion 
we stress that there is no parallel between the two formulas Wr'~'^ (T) and W^r^(r,ro): on 
the one hand Wr'~'^ (F) yields the correct value of the writhe of a curve as soon as the curve 
does not self-intersect and on the other hand Wr^{T, Fq) often yields an incorrect value for 
the writhe of a curve even if this curve has no antipodal point. 

All what is exposed here for the curves in Fig.[7]is directly applicable to the configurations 
of Fig. [2] and shows that the writhe (resp. the link) of the configurations with Z < Zb in 
Fig. His not given by Eq. Q (resp. Eq. ©). 

Impossibility to use the 'stadium' shaped curve as reference curve 

Fuller's theorem [T^] states that as soon as there is a continuous deformation between F and 
Fq which is free of antipodal points and self-crossings (we called 'good' such a deformation), 
VFr^(F,Fo) yields the correct value for the writhe of the closed curve F. In Fig. [7] the 
deformation from A = OtoA = lisa 'bad' deformation as it contains one antipodal event. 
We then conclude that PFr^(Fi, Fq) yields an incorrect value. But one might say that another 
deformation between Fi and Fq, with neither self-crossing nor antipodal event, could exist. 
We show here that this is not the case. Statement: no 'good' deformation between Fi and 
Fq of Fig. [7] exists. Proof, assume a 'good' deformation exists. Then we can build a closed 
circuit by joining this proposed deformation and the (reversed) actual deformation A G [0, 1] 
in Fig. [71 One would then have a closed circuit with one antipodal point and no self-crossing, 
in contradiction with the result established in [22] . Hence no such 'good' deformation exists. □ 
This means that as soon as one 'bad' deformation exists between a reference curve Fq and a 
curve F the formula M^r^(F,Fo) does not yield the correct result and no 'good' deformation 
can exist. A consequence of this is that rotating the reference curve of Appendix A (Fig. [6l 
top) to introduce a new, and 'good', reference curve (e.g. Fig. bottom) is hopeless. To 
show this we first introduce a 'bad' deformation between the rotated reference curve (e.g. 
Fig. [HI bottom) and the supercoiled plasmid of Fig. O we untangle the plasmid using a self- 
crossing, as in Fig. [7]-C, to obtain a 'stadium' shaped curve which we subsequently rotate to 
the (rotated) reference curve. The existence of this 'bad' deformation means that no 'good' 
deformation can exist. 
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